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Computing the entanglement of formation of a bipartite state is generally difficult, 
but special symmetries of a state can simplify the problem. For instance, this allows 
. one to determine the entanglement of formation of Werner states and isotropic states. 

We consider a slightly more general class of states, rotationally symmetric states, also 
known as SU(2)-invariant states. These states are invariant under global rotations of 
both subsystems, and one can examine entanglement in cases where the subsystems 
have different dimensions. We derive an analytic expression for the entanglement of 
formation of rotationally symmetric states of a spin-j particle and a spin-| particle. 
, We also give expressions for the I-concurrence, I-tangle, and convex-roof-extended 

£> \ negativity. 
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I. INTRODUCTION 



An entangled pure state is one in which complete knowledge about the overall state is 
incomplete with regard to the subsystems Such states are strongly correlated in the 
sense that the correlations cannot be reproduced by a local hidden variable theory 0, E| . 
Besides being integral ingredients in no-go theorems regarding the interpretation of quan- 
tum mechanics, entangled states open new avenues in quantum engineering and information 
processing (Ula]. Using entangled states, we can teleport quantum information, we can in 
principle build quantum computers that provide exponential speedup over classical comput- 
ers, and we can also prove the security of assorted cryptographic protocols p. We can do 
all this and more without a complete understanding of entanglement, because the tasks rely 
only on pure states. Yet one often must deal with mixed states that arise from physical 
processes such as noise. This makes it important to have a measure of entanglement for 
mixed states so that potential applications can be explored and evaluated. 

A commonly used entanglement measure for a pure state of two systems, A and B, 
is the entropy of the marginal density operator pa (or p B ), 

E(ip) = S(p A ) = -tr(p A logp A ) = -ti(p B \ogp B ) . (1.1) 

The importance of this measure comes chiefly from the fact that it gives the rate at which 
copies of a pure state can be converted, by using only local operations and classical com- 
munication (LOCCs), into copies of maximally entangled states and vice versa 0,111. The 



entanglement of formation 0, 

E F (p) = min l^PjEty 

\ 3 



is the so-called convex-roof extension of the pure-state measure (jl.l|) to bipartite mixed states 
p. The need to search over all ensemble decompositions of p generally makes it impossible 
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to calculate the entanglement of formation, unless some efficient method to short-circuit a 
complete search can be found. The entanglement of formation provides an upper bound 
on the rate at which the maximally entangled states must be supplied to create copies of 
P 0> 11 an d the rate at which one can distill maximally entangled states from p. It is also 



a measure for which analytic results are known. The entanglement of formation is known for 
arbitrary bipartite qubit states [13], isotropic states [3] and Werner states in arbitrary 
dimensions, and symmetric gaussian states in infinite dimensions [IsT ]. 

Many other mixed-state entanglement measures have been defined: entanglement of dis- 
tillation 0, negativity EB, relative entropy [l7j . robustness of entanglement [IH, I- 



concurrence |19I |. I-tangle |20], the geometric measure 2l|, and others. Each is useful in 
particular physical contexts, but the different measures produce different orderings on mixed 
states even when they agree on pure states 2^1 . Except for negativity, these measures are 



also difficult to calculate for an arbitrary mixed state. Here we focus on the entanglement 
of formation, but provide results for other measures since our results on entanglement of 
formation can easily be extended to these other measures. 

Terhal and Vollbrecht [13] showed how to use the symmetry of the isotropic states to 
find the entanglement of those states. Vollbrecht and Werner [l4( elucidated and extended 
the results of that paper; in particular, they formulated a general technique that leverages 
symmetries to simplify the calculation of convex-roof measures for symmetric states. We 
apply their method to a bipartite system consisting of a spin-j particle and a spin-| particle, 
i.e., a [(2j + 1) x 2]-dimensional Hilbert space, where the states we consider are those that 
are invariant under global rotations of the two particles. Such states are also known as 
SU(2)-invariant states jl^]. They are functions of a single parameter p, 



P(P) = ^~^ n j+ 1 /2 + f- n i-l/2 > ( L3 ) 



where the operators 



i±i/2 



n i± i/ 2 = 2^ |j±|,m)(j±i,m|, (1.4) 

m=-(i±l/2) 

are the projectors onto the subspaces of total angular momentum j ± |. These states can 
appear when one loses information regarding the Cartesian reference frame of the two sys- 
tems |24|. They can also arise from bipartite splits of a rotationally symmetric chain of 
spin- 1 particles into a single qubit and the rest. Such states can also arise as the multi- 
photon states that are generated by parametric down-conversion and then undergo photon 
losses jiil ]. 

The rotationally symmetric state ()1.3|) is known to be separable if and only if 2j| 2(| 



p < 2j . (1.5) 
2 j + 1 y } 

In this paper we show that the entanglement of formation of the SU(2)-invariant state (jl.3|) 
is given by 

pe [0,2j/(2j + l)], 

M ' H = l » / ^ I (v5-«)| Pe [2,7(2, + (L6) 
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Here H(x) = — xlogx — (1 — x) log(l — x) is the binary entropy; we use the natural logarithm 
in calculating actual values of the entanglement of formation and when taking derivatives of 
entropic expressions. The formula (jl.6|) is noteworthy in that it provides the first example 
of an entanglement of formation for subsystems having different dimensions. 

In Sec.|nJ we briefly review the Terhal-Vollbrecht- Werner program for determining the en- 
tanglement of symmetric states. Section ITTTI considers rotationally invariant states and their 
properties. In Sec. IIV| we calculate the entanglement of formation of rotationally invariant 
states of a spin-j particle and a spin-| particle. Section |V| contains analytic expressions for 
other mixed-state entanglement measures. Finally, in Sec. IVII we address possible extensions 
of our results. 



II. CONVEX-ROOF ENTANGLEMENT MEASURES UNDER SYMMETRY 

In this section we summarize the Terhal-Vollbrecht-Werner procedure for determining; 
the entanglement of symmetric states. 

We begin with a few definitions. Let K be a compact convex set (e.g., density matrices), 
which is itself a subset of a finite-dimensional vector space V, and let M C K be an arbitrary 
subset of K (e.g., pure states). Let / : M — > RU{+oo} be a real-valued function on M (e.g., 
an entanglement measure on pure states). The convex roof of /, cof : K — > R U {+00}, is 
a function on the entire set K, defined by 



cof(x) = min > A,f(s,) 



$>^ = x, A,>0, = ! . (2.1) 



The convex roof is computed by minimizing the average value of / over all possible ways 
of writing an element x G K as a convex combination of elements in M. In the case of 
density operators and pure states, such a convex combination is called a pure-state ensemble 
decomposition, so we refer generally to a decomposition of x. If x has no decomposition 
in terms of elements of M, cof(x) is infinite. The convex roof cof(x) is the largest convex 
function g on K such that g(s) < f(s) for s G M; if M is a set of extreme points of K, as 
in the case of pure states and density operators, co/(s) = f(s) for s G M. In this notation, 
the entanglement of formation (|1.2j) is 

E F (p) = co E(p) . (2.2) 

Now suppose there exists a symmetry group Q that acts on V through its matrix repre- 
sentations olq- Assume that olqK C K and also that agM C M, with /(«g(s)) = f(s) for 
all s G M. In the case of entanglement, the linear transformations «g are tensor products 
of local unitary representations of a symmetry group; as unitary transformations, they take 
pure states to pure states, and as tensor products of local unitary transformations on the 
two subsystems, they preserve the entanglement measure. Define a projection P : K — > K 
by averaging uniformly over the group, 

P(x) = [ dGa G (x) . (2.3) 



The projection P is often called twirling. We denote the range of P by PK. Elements of PK 
are precisely those elements of K that are invariant under the projection. More generally, 
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PK consists of all elements of K that are invariant under the group operations and thus are 
known as group-invariant elements. It is clear that PK is a convex set. 

The problem of finding cof(x) for x G PK is simplified, because the problem can be 
divided into two parts. One first determines the function e : PK ->1U {+00} defined by 

e(x) = min (f(s) I P(s) = x) . (2.4) 

{seM} y 1 

Thus instead of minimizing over all possible decompositions, one first minimizes over a 
subset of the pure states. Then the convex roof of / on PK is the convex hull of e : 

cof(x) = coe(x) , x G PK . (2.5) 

The convex hull of e is the largest convex function on PK that nowhere exceeds e. We use 
the same notation for convex hull and convex roof, because the convex hull of a function / 
is the convex roof for the special case M = K, i.e., in this case, 



coe(x)= mm > A,elx,- 



J2 = ^ > 0, Yl x > = 1 ) 

3 3 j 



(2.6) 



For highly symmetric states, the more difficult of the two steps is determining e(x), since 
finding co e(x) is usually straightforward or even unnecessary because e(x) is already convex. 

Demonstrating the reduction of the minimization (j2.1|) to Eqs. f!2.4|) and ()2.5|) is suffi- 
ciently simple that we include it here for completeness. First, let x G PK have the opti- 
mal decomposition x = J2j X j s j relative to Eq. (|2.1j) . i.e., cof(x) = J2j X jf( s j)- Defining 
Xj = P(sj), we have e(xj) < f(sj), x = P(x) = ^ • XjXj, and 

co/(x) = X jf( s i) ^ E Vfo) ^ co < x ) ■ ( 2J ) 
3 j 

Second, let x G PK have the optimal decomposition x = J2j X j x j relative to Eq. (|2.6|) . i.e., 

coe(x) = Y2j x j e ( x j)i ari d let Sj achieve the minimum in Eq. ()2.4j) for Xj, i.e., Xj = P(sj) 
and e(xj) = f{sj). Then 

x = ^2\j / dGa G {sj) (2.8) 
3 ^ 

is a decomposition of 2, and 

co/(x) <J2 X j[ dGf(a G { Sj )) =J2 X jf(s 3 ) = J2 X J< x i) = co < x ) ■ (2-9) 

Together, Eqs. (j2Hj) and (Q give Eq. (I23J) . 

This method has been used to compute the entan glem ent of the Werner states h|| , which 
are invariant under U ® U, and the isotropic states |l3j], which are invariant under U <8> C/* • 
One can generalize to a larger class of states by considering a subgroup of either U ® C/ or 
U ®U* . In the next section we examine groups of the form i? ® i£, where i? is a rotation. 
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III. SU(2)-INVARIANT STATES 

Consider two particles, one of spin ji and the other of spin j 2 . We consider states that 
are symmetric under global rotations. The symmetry group is 

where D^(R) = exp(—i8J ■ n) is the spin-j representation of the rotation R e 5*0(3) [or 
equivalently in SU(2)]. Here J is the angular- momentum vector, with components J x , J y , 
and J z . The twirling operator is then 

P n (p) = J dfi{R) D {jl) (R) <g> D (j2 \R)pD {jl) (Ry ® D {n) {R) ] , (3.1) 

where fi{R) is the group-invariant measure for the rotation group. The twirling operation 
describes a process where two parties use classical communication to select a random rotation 
R that each party implements locally. This makes twirling a LOCC operation, which implies 
that entanglement does not increase under twirling. 

The states that are invariant under the twirling operation are those that are convex 
combinations of the states associated with the projectors onto the irreducible subspaces of 
total angular momentum J = \ji — j 2 | , \ ji — J2I + 1, • • • , ji + ji- These projectors, 

j 

n J= \J,m)(J,m\ , (3.2) 

m=—J 

have associated normalized states Ilj/(2J + 1). Thus the 7?.-invariant states are 

ji+h 

*>>= £ wU*^ (3-3) 

J=\3l~fe\ 

where pj = tr(p(p)Ilj) > and ^2pj = I- Any state a twirls to a 7£-invariant state, i.e., 
Pr((j) = p(p), where pj = tr(P^(cr)nj) = tr(o"P-^(rij)). Since ITj is 1Z invariant, we find 
that pj = tr(crllj) is the overlap of a with the subspace of total angular momentum J. 
Notice also that 

P n (\J,m)(J,m\) (3.4) 

for any value of m. 

We note that, for the ^-invariant state ()3.3j) . the positive partial transpose condition is 
necessary and sufficient for separability when j\ — \ and j'2 is arbitrary and when j\ — 1 
and j2 is an integer 0, Q H . 

The general problem of finding the entanglement of formation of rotationally invariant 
states is first to determine the function 

e(p) = minto) P^l^^l) = p(pj) = min(^) | (m^) = Pj) (3-5) 



{|^>}V 



{!'/■>} 



and then to find its convex hull on the convex set of probabilities p. 

It is interesting to note that although states invariant under U ®U and U ®U* are quite 
different, this is not the case if one compares states invariant under i? eg) R and R® R* . 
The reason is that conjugation of a representation is equivalent to rotating by tt about 
the y axis, i.e., D^\R)* = e lnJy D^\R)e~ tnJv . Thus the states that are invariant under 
1Z = {D^ 1 \R) ® D^ 2 \R)*} are obtained from the rotationally invariant states by rotating 
one of the two systems by 7r about the y axis. 
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IV. ENTANGLEMENT OF FORMATION FOR SU(2)-INVARIANT STATES OF 

SPIN-j AND SPIN- 1 PARTICLES 

A. General considerations 

To obtain insight into the entanglement of rotationally symmetric states, we examine the 
simplest case, a particle of arbitrary spin j and a particle of spin-|. The eigenstates of total 
angular momentum are well known, given by the Clebsch- Gordon coefficients |3CI| . and can 
be written as 



\j ± h m) = ± J '-^f \J, ™-l)®\H) + J ^±f^ \j, m+\)® li ~\) ■ (4.1) 



It is straightforward to calculate the entanglement of these states: 

£(u±i , ro) )^(i±i^)^(i-J^). ( , 2) 

Eigenstates with identical values of \m\ have the same entanglement. In the j + \ subspace, 
the minimum entanglement is 0, achieved when \m\ — j + ~ and in the j — \ subspace, the 
minimum entanglement is if (1/(2 j + 1)), achieved when \m\ = j — |. 

The ^-invariant states are convex combinations of the states associated with the projec- 
tors ()1.4|) onto the subspaces of total angular momentum j + ~ and j — \, 

Pip) = 2{] + l) Y[3+1/2 + 2j Uj ~ 1/2 ' ^ 4 ' 3 ^ 

where p = tr(p(p)nj_i/2) . Any state a twirls to P-^cx) = p(p), where p = tr(crITj_i/2) is the 
overlap of a with the j — \ subspace. 

The first reason this problem is simpler than the general case of two arbitrary spins is 
that the rotationally invariant states are specified by the single parameter p. The problem 
of finding the entanglement of formation reduces to determining a function of this one 
parameter, 

e(p) = min( E(ip) \ (^-1/2^) = p) , (4.4) 
{lw 

and then finding its convex hull, coe(p). 

B. Determining e{p) 
To find e(p), we begin by looking at a couple of example states. The first, 



10) = |j',J>®(VT^|i,i) + v^l5>-5>) 

= v^|j + ij + |) + y^|j + |,J-|) + y^|j-|,J-|), (4.5) 

is a product state, thus having no entanglement, i.e., E(<p) = 0. This state has overlap 
p = (</>|IL,-_i/2|</>) = 2ju/(2j + 1). As v ranges from to 1, p varies from to 2j/(2j + 1), 
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which shows that e(p) = for p G [0, 2j / (2j + 1)]. Thus in determining e(p), we can restrict 
our attention to p G [2j/(2j + 1), 1]. 
The second state, 

Ix) = -Vi"Ui- 1 )®ll>l> + y/r=7i\jJ)®\h-k) 

' f (-V / 2^+ V / W)b' + |,J - I) + >/2j(l-/x))|j - |, j - §) 



(4.6) 



is entangled, with J5(x) — H(fj), and has overlap 

p = (xlKj-iMx) = 2j~+i (v^+ VW^)) 2 = P„ ■ (4-7) 

As /x increases from 0, p increases monotonically from a value of 2j / (2j + 1) at /x = to a 
maximum value of 1 at /x = 1 / (2j + 1); for larger values of /x, p decreases monotonically to 
a value of l/(2j + 1) at /x = 1. 

Inverting to find /x as a function of p and using the branch that gives the smaller values 
of /x, we find that for p G [2j/(2j + 1), 1], 

V = ^7 (VP " \/2j(1-P)) 2 = /imin(p) • (4.8) 

The reason for the functional notation /x m i n (p) becomes clear below. The function /x m i n (p) 
increases monotonically from a value of /x = at p = 2j/(2j + 1) to /x = 1 / (2j + 1) < 1/2 
at p = 1. The upshot is that the state (|4.fij) tells us that for p G [2j/ (2j + 1), 1], 



e(p) < ^(/i min (p)) = ^(^t(vp- >/2j(i -p)) 2 ) < ^(aj^T 



(4.9) 



It turns out that |x) achieves the minimum value of E(ip) in Eq. f)4.4j) and thus that e(p) 
is actually given by the expression in the middle of Eq. (|4.9jl . To proceed with the proof of 
this, however, the only information we need from \x) is the final inequality in Eq. (|4.9j) . i.e., 
e(p)<F(l/(2j + l)). 

The next step is to characterize the set of pure states over which we must minimize in 
Eq. ([4.4)1 . Any pure state of a (2j + 1) x 2 system can be written in terms of a two-term 
Schmidt decomposition, 

W) = v^lei) ® l/i) + v / W|e 2 ) ® |/ 2 )) • (4.10) 

which has entanglement -E(V') = H(p). The second reason this problem is simpler than the 
general case is that the pure states that twirl to invariant states are specified by the single 
parameter /x. The final inequality in Eq. ()4.9|) implies that the only values of /x we need to 
consider are /x G (0, 1 / ( 2 j + 1)], for which H(fi) is monotonically increasing. 

We can introduce two unitary operators, V on the spin-j particle and W on the spin-| 
particle, which transform bases of our choice to the Schmidt bases. In particular, we can 
write j^) as 

|^) = V ® W\ X ) = V® W(-Jji - 1) <g> |i, \) + VT=]i ® li -i)) • (4.11) 
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Any unitary on the spin-| system is equivalent (up to an irrelevant phase) to a rotation R, 
i.e., W = D^' 2 \R). Rotating both particles by the inverse of R, we obtain the state 

\MU)) = u® i\x) = u® /(-v^ - 1) ® \h \) + V^\j,j) ® \h 4» > ( 4 - 12 ) 

where C/ = D&')(fl)tv\ 

The significance of this move is that \ip^{U)) has the same entanglement and the same 
overlap with the j — | subspace as does Thus in doing the minimization ([4.4)1 . we only 
need to consider the states |?/> M (f/)), i.e., 

e(p) = mm (H(p)\p,(U)=p) , (4.13) 

where 

P,(U) = (%(U)\U^ 1/2 \^(U)) = ( X \rf®IIl 3 _ l/2 U®I\x) . (4.14) 

Moreover, since H(p) is a monotonically increasing function of \i in the interval of interest, 
we can simply minimize fi over all unitaries U, 

jl(p) =nun(/i | p^U) =p) , (4.15) 

and plug the result into the binary entropy to give 

e(p) = H{p) . (4.16) 

The final reason that the problem of finding the entanglement of formation in this case is 
doable is that the problem can be reduced to the minimization ()4.15|) over a single unitary, 
that being a unitary acting on the spin-j particle. 

Our strategy now is to show that p^{U) < PfJ ,(I) = p^, where p^ is the function (|4.7|) . 
This result, which we will prove below, allows us to immediately determine p.(p). Since p M is 
monotonically increasing on the relevant interval (0, 1 / ( 2 j + 1)], we can find a // < \i such 
that p^(U) = Pn'(I) for any p and U, implying that the minimum is always achieved by 
U — I. This means that p is obtained by inverting p^ = p, and thus 

HP) = /^min(p) , (4.17) 

as given by Eq. ()4.8jl . As promised, e(p) on the interval [2j / (2j + 1), 1] is given by the middle 
expression in Eq. ()4.9|) . The result is that 

(0, pe[0,2j/(2j + l)], 

€(p) = \»(^(v-p-VM^7)) 2 ) , P e m2j + i),i]. (418) 

We now must show that p^U) < p^. We begin by noting that 

MU) = \(x\tf®i\j-\M\ 2 

m=-j+l/2 

= ^T^TT S | V 7 ! 1 - V)r m x m e iam + ^JIF^y m e if3m , (4.19) 

m=-j+l/2 
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where we define two rows of matrix elements of W, 

{j,j\tf\j,m + §) = x m e iam , x m > , m = -j - \, . . . ,j - \ , (4.20) 

- l|E% m _ I) = ym e^ , y m > , m = -j + |, . . ., j + | , (4.21) 

and where 

r m =j + l + m>0. (4.22) 

Only two rows of the unitary matrix are involved in the overlap (J4.19j) because one particle 
is a qubit. Now we can write 

fV(^) - 2 ■ + i _ t^VrnXm + ^nr„ m y m ^ = - , - , (4.23) 

m=— j+1/2 

with equality when all the phase factors are equal to 1. For the remainder of the proof, we 
omit the range of the sum since it is always that of Eq. ()4.23|) . i.e., m — — j + |, . . . ,j — 
Notice that on this range, r m and r_ m are strictly positive. 

The problem now is one of maximizing F(x, y) subject to the constraints on the two 
relevant rows of W. Normalization of these rows gives the constraints 

Yl = 1 - X -j-l/2 and Yl ^ra = 1 ~ • ( 4 - 2 4) 

The two rows must also be orthogonal, but we ignore this constraint on the grounds that 
doing so can only lead to a bigger maximum, which is still an upper bound for p^{U). With 
this constraint ignored, it is clear from the form of Eq. ()4.23j) that the maximum is achieved 
when X-j-1/2 = = l/j+i/2- Thus we introduce two Lagrange multipliers and maximize the 
function 

G(x, y) = F(x, y) + Ai [J2 *l ~ l) + A 2 y 2 m - l) - (4.25) 

It turns out that solution to this maximization satisfies the orthogonality constraint. 
Setting the first derivatives equal to zero yields the following equations, 

[(1 - fi)r m + Xi]x m + a/(1 - fi)fir m r_ m y m = , 

a/(1 - fi)fir m r_ m x m + (jir-m + A 2 )y m = , (4.26) 

which hold for all m. Since x m and y m are both nonnegative, we must have Ai < — (1— \i)r m < 
and A2 < —fir^ m < 0. Moreover, if the solution is not to be x m = y m = 0, the determinant 
of the matrix of coefficients must vanish, which gives 

A X A 2 + (1 - i2)r m X 2 + /xr_ m A a = . (4.27) 

Now suppose Eq. ()4.27|) holds for two or more values of m. By considering any pair of m 
values for which Eq. ()4.27|) holds, it is easy to show that the Lagrange multipliers are given 
by Ai = A2 = or by 

Ax = -(1 - n)(r m + r_ m ) = -(1 - /i)(2j + 1) , (4.28) 
A 2 = -nir m + r_ m ) = -n(2j + 1) . (4.29) 
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The former case is ruled out by the requirement that Ai and A2 be negative, so we need only 
consider Eqs. (14.28(1 and ()4.29j) . which imply 



X m = \ Z \ — Vra ■ (4.30) 

V 1 - fx V r_ m 

Now we find that 

E^T^Er^r^EA* 1 . < 4 - 31 > 

with equality if and only if the only nonzero y m is m = j — | and p = 1 / (2j + 1). We 
conclude that it is impossible to have nonzero solutions for x m and y m for more than one 
value of m. 

The result is that the only solutions of the derivative equations ()4.26() have just one 
value of m for which x m and y m are nonzero, and for that value, the constraints imply that 
%m — Vm — I; giving 2j extrema of F(x,y). For all these extrema, the two rows of W are 
orthogonal, as required by the constraint we neglected. The value of F at the extreme points 
is a function of m, 

F(x, y) = F m = (y/(l - /i)r m + ^pTI^) 2 . (4.32) 

For < p < l/(2j + 1), F m is monotonically increasing for m e [—j + |, j — |]. Hence the 
maximum occurs at m = j — |. Pulling all this together, we have 

P,(U) < yJjFj-i/* = 2j^T (v^+ V2JM) 2 = P, ■ (4-33) 
This establishes the result we needed above to complete our determination of e(p). 



C. Determining co e(p) 

We now have to find coe(p). What we show is that e(p) is convex, so coe(p) = e(p). 
The second derivative of e(p) is given by 

W , - p) )^e» W = ^ log ( ^; 2 ^! ? p) ) - ■ (4.34) 

We plot e"(p) for j — |, 1, and 3 in Fig. HI and we see that the second derivative is positive 
for these values of j. This turns out to be true for all j, as we see in the following way. It is 
straightforward to show that the right-hand side of Eq. (|4.34j) is a nonincreasing function in 
the region of interest by showing that the derivative is bounded above by zero, and therefore 
its minimum value occurs at p — 1. We thus get the result 

(p(l - p)f /2 e"{p) > 2( ^ 1} log2j , (4.35) 

and so the second derivative of e(p) is always positive for j > 1/2, which establishes that 
e(p) is convex. 

We have established that the entanglement of formation is given by Epipip)) = coe(p) = 
e(p), thus verifying Eq. (|1.6(1 . We plot the Ep(p{pj) for j = |, 1, and 3 in Fig. 
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log (6" (p) ) 




! 6 0.1 0.8 ! 9 1 P 

FIG. 1: loge"(p) for j = | (solid), j = 1 (long-dashed), and j = 3 (short-dashed). We plot the log 
because e" diverges rapidly at p = 1 for all j > 1. Since the log is positive, the second derivative is 
bigger than 1, and therefore e{p) is convex. Equation (|4.35f) shows that this holds for all j > 1/2. 

E(p(p))=e(p) 




FIG. 2: The entanglement of formation in units of nats (we use the natural log in the binary 
entropy) for j = \ (solid), j = 1 (long-dashed), and j = 3 (short-dashed). The rotationally 
invariant states become less entangled as j increases. 



V. CONVEX-ROOF MEASURES 

Many entanglement measures can be constructed using the the convex-roof extension. 
Given any pure-state entanglement measure X{ip) = f(pA) that is (i) invariant under local 
unitaries and (ii) a concave function of the marginal density matrix p^, coX(p) is an entan- 
glement monotone [31]. The procedure used here to find the entanglement of formation of 
rotationally invariant states applies to any other convex-roof extension: first determine 

£(p) = mm(X(^) | (Vln^lV) =p) , (5.1) 
{\i>)} 

the minimum value of X over all pure states that project to p(p), and then compute the 
convex hull of £(p). The state |%) of Eq. ()4.6|) . which minimizes e(p), also minimizes £Qo), 
because different pure-state entanglement measures give the same ordering of pure states. 
In terms of the quantities introduced in our consideration of the entanglement of formation, 
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this is the statement that for the state ([4.10)1 . /(//) = X{jp) is an increasing function of \i 
for < \i < 1/(2 j + 1), thus giving = /(p m m(p)), in analogy to Eq. ()4.16jh 

In this section we apply this technique to find the concurrence, the tangle, and the 
negativity of rotationally invariant states of a spin-j particle and a spin-| particle. 



I- Concurrence 



The generalized concurrence of a joint pure state of systems A and B measures 
the purity of the marginal states. 



Cty) = V 2 [1 - tr(pi)] = V 2 t 1 - tr (p|)l • ( 5 - 2 ) 

The concurrence of the state (jCTIJ) is C{ip) = 2^(1 - //). The convex-roof extension of 
the concurrence, C(p) = coC(p), is called the I- concurrence. 

As discussed above, our results for the entanglement of formation determine the function 

c(p) = miri(C(V>) | (V>|n i _i/ !! |V') =P) 
llw 

2 ^/pmin(l pmin) 

= — (V27(2p - 1) - (2j - 1)^(1 - p) ) • (5.3) 

The function c(p) is clearly convex because it is the sum of convex functions. When both 
particles have spin-i, c(p) is linear, c(p) — 2p — 1 for 1/2 < p < 1, and it becomes more 
convex as the spin j of the first particle increases. Thus the I-concurrence of a rotationally 
symmetric state is 

- \c(p) , p £[2,7(2; + !),!]. (5 ' 4) 
The I-concurrence increases from at p = 2j/ (2j + 1) to 2^/2] / (2j + 1) at p = 1. 



B. Tangle 

The tangle is the convex- roof extension of the squared concurrence [20] : 

r(p) = co C 2 (p). (5.5) 

Since 

min(C 2 W | (mj-i/M) =P) = c 2 (p) (5.6) 
llw 

and the convexity of c(p) implies the convexity of c 2 (p), the tangle of rotationally symmetric 
states is given by 

[P[P)) ~ \c\p) , pe [2j/(2j + !),!]. (5 ' 7) 
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C. Convex-Roof-Extended Negativity 

The negativity of a joint density operator p is proportional to the sum of the negative 
eigenvalues of the partial transpose of p: 

Af(p) = -2J2^ = \\P TB \\-1- (5-8) 

Here the p/s are the eigenvalues of the partial transpose p TB of p with respect to system 
B, and ||p Ts || is the sum of the absolute values of these eigenvalues. The negativity is an 
entanglement monotone and it is straightforward to evaluate. The negativity does not, 
however, identify all entangled states. There are states, the bound entangled states, that 
have a positive partial transpose and thus zero negativity, yet are still entangled. Notice 
that there are no such states in the one-parameter family of rotationally symmetric states. 
To get around the inability of the negativity to identify all entangled states, the convex- 
roof-extended negativity, coAf, has been proposed as a mixed-state entanglement measure. 
This comes with a considerably increased difficulty in computation, of the sort associated 
with any convex-roof entanglement measure, but the convex-roof-extended negativity has 
been evaluated for the isotropic states and Werner states [32| . 

For p e [2j/ (2j + 1), 1], the partial transpose of the rotationally symmetric state p{jp) of 
Eq. (jl.3|) has two eigenvalues 

p_ = — , (5.9) 

^ 2j + l 2j' 1 ; 

the first of which is 2(j + l)-fold degenerate, and the second of which is 2j'-fold degenerate 
and is negative for p e (2j/(2j + 1), 1] |23j. Thus the negativity of p(p) is 



max 



0.21/; 2j 



2j + l 



(5.10) 



To find the convex- roof-extended negativity, we note that for the states ()4.10|) . J\f(ip) = 
2^/ p(l — p) = C(ijj). This means that for rotationally symmetric states, the convex-roof- 
extended negativity is actually identical to the I-concurrence (15. 4|) . 



CO 



Af(p(p)) = co C(p(p)) = C(p(p)) . (5.11) 



VI. SUMMARY 

We have computed the entanglement of formation of rotationally symmetric states for a 
[(2j + 1) x 2] -dimensional system. Three features of the problem simplified the endeavor: the 
rotationally invariant states are specified by a single parameter, pure states that are twirled 
to an invariant state are determined by a single constraint, and only one local unitary is 
needed to relate the pure states to a fiducial state. The first makes finding the convex 
hull in Eq. (j2.5|) easier, and the other two allow one to efficiently characterize the subset 
of pure states over which one minimizes in Eq. (j2.4j) . These elements were also vital in 
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all other applications of the Terhal-Vollbrecht- Werner procedure to date, i.e., generalized 
Werner states and isotropic states. 

The rotationally symmetric states of a spin-j particle and a spin-| particle are similar 
in some ways to the states of a two-qubit system: positive partial transpose is necessary 
and sufficient for separability, and the pure states R ® R\x) i n the optimal decomposition 
all have the same entanglement. The entanglement in these states vanishes as one spin 
becomes more classical (j — > oo). We have also determined the I-concurrence, tangle, and 
the convex-roof negativity of the rotationally symmetric states. These all display the same 
features as the entanglement of formation. This is quite different from what occurs for the 
isotropic states in high dimensions, where the concurrence is linear [33;], the tangle has the 
same behavior as the entanglement of formation [33|, and the convex- roof negativity is the 
same as the negativity of the isotropic state • 

Since we have determined an optimal decomposition, we can also find the entanglement 
of formation of states that are not rotationally symmetric. This procedure for extending 
optimal decompositions to other states was outlined by Vollbrecht and Werner [l4|]: arbi- 
trary convex combinations of the pure states R <S> R\x) have the same entanglement as the 
rotationally symmetric states, which are uniform convex combinations of these states. This 
extension procedure does not cover the entire space; the states for which it does apply and 
their properties are currently being investigated. 

We have examined the simplest case of rotationally invariant states, and the next step, 
which we are pursuing, is to investigate two spin-1 particles, where the rotationally sym- 
metric states constitute a two-parameter family. For two spin-1 particles, it is known that 
the positive partial transpose condition is necessary and sufficient for separability, and the 
entanglement of a large fraction of the rotationally symmetric states can be obtained from 
extending the decompositions for generalized Werner states and isotropic states 0]. How- 
ever the calculation for the rest of the states poses some difficulties, because the important 
features that facilitated the calculations in this paper are now absent. The states are func- 
tions of two variables, there are two constraints on the pure states, instead of one, and it 
becomes difficult to characterize efficiently the states that are twirled to an invariant state. 
In particular, the technique used to remove the unitary on one subsystem, thus obtaining 
Eq. (|4.12|) . no longer applies. The entire problem gets progressively more difficult as the 
spins of the two particles increase. 
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